AN INVISCID DYADIC MODEL OF TURBULENCE: THE 
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Abstract. Properties of an infinite system of nonlinearly coupled or- 
dinary differential equations are discussed. This system models some 
properties present in the equations of motion for an inviscid fluid such 
as the skew symmetry and the 3-dimensional scaling of the quadratic 
nonlinearity. It is proved that the system with forcing has a unique equi- 
librium and that every solution blows up in finite time in _ff 5 ^ 6 -norm. 
Onsager's conjecture is confirmed for the model system. 



1. Introduction 

One of the outstanding questions in fluid dynamics is existence, unique- 
ness and regularity of solutions to the Cauchy problem for the three-dimensional 
Euler equations: 

(1.1) ^ = -(«-v)«-v Pl 

V-« = 0. 

Important features of the nonlinear term (u ■ V)u are its bilinearity and 
skew-symmetry. The second property implies conservation of energy for suf- 
ficiently regular solution of the Euler equation (and decay of energy in the 
context of the Navier-Stokes equations). In the past few decades discretized 
"toy" models that preserve the energy properties of the fluid equations have 
been proposed and studied both by mathematicians and physicists. These 
models belong to a general class of "shell" models which simulate the en- 
ergy cascade in turbulent flow. In all of these models the nonlinearity of 
the Euler equations is much simplified by considering only local neigbouring 
interactions between certain scales. However simplifications vary and as a 
consequence the models differ in the number of conserved quantities and in 
the presence of a certain "monotonicity" property that we will discuss later. 
Among the first example of such discretized models is the one introduced by 
Gledzer |llj in 1970 which was later generalized by Ohkitani and Yamada 
|17j and is now known as the GOY model. A survey of mathematical devel- 
opments in connection with shell models can be found in the recent book of 
Bohr et al [T]. 
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In this article we analyze one of these discretized models, namely 
^ = 2*G- 1 >a?_ 1 - 2§ + fr, i>o 

(1.2) , 

— = -a ai + / . 

The motivation for this model, and specifically for the scaling appropriate 
for 3 dimensions is given in [10] and [12] in the case —00 < j < 00 and 
/ = 0. In [10] a shell model is presented, following the Fourier space anal- 
ysis of Dinaburg and Sinai [8], for the 3-dimensional incompressible Euler 
equations (jl.lj) and it is shown that for very specific initial data the vec- 
tor model reduces to the scalar system analogous to (|1.2p without forcing. 
Furthermore, the wavelet decomposition utilised in [12] and [13] motivates a 
discretized model for the 3-dimensional Euler equations analogous to (|1.2p 
without forcing. 

A derivation of the system (jl.2p is presented in Section [2j The coefficient 
a 2 (t) is the total energy in the frequency space shell 2^ < \k\ < 2 3+1 . In this 
context I 2 and H s , respectively energy and Sobolev norms, are defined as: 



1* 



y/ 2 (00 \w 

=0 / \j=0 



Local in time existence of solutions in H 5 / 2 was obtained in [10] using the 
fixed point techniques that produce local existence in H 5 ^ 2 for mild solu- 
tions of actual 3 dimensional Euler equations (11. ip (see Cannone [2]). Finite 
time blow-up of (|1.2j) . without forcing, in an appropriate H s was initially 
proved in [12] by exploiting conservation of energy and a monotonicity prop- 
erty present in the quadratic term Variants of the model have been 
studied in [T3], [20], and [21]. In particular, it was observed in [20J that 
(jl.2p with 2^l 2 replaced by 2- 7 could be derived from a discretized version 
of the Fourier transform of the 1-dimensional Burger equations under the 
(unrealistic) restriction that there is no spreading of the support. 

Variants of the model including dissipation have been studied [3], [7], 
[12] . [13] and finite time blow up is proved for certain "small" amounts 
of dissipation. More precisely, in [3] a finite time blow up was proved in 
the cases where the dissipation degree is such that the model enjoys the 
same estimate on the nonlinear term as the Navier-Stokes equations in the 
dimensions larger than 4. 

Even though a global existence proof from [3] applied to (11.2(1 implies that 
for any initial data in I 2 there exists a global in time classical solution to 
(|1.2p . not all the classical solutions satisfy the energy balance equation. For 
instance, regular solutions, i.e., solutions with bounded i? 5 / 6 -norm, satisfy 
the energy equality, but the unique fixed point does not. 

In this present paper we study the system (|1.2p with forcing as a model 
for 3-dimensional turbulence. The main results are the following: 
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(a) For any force / = (/o,0, ...), /o > 0, there exists a unique fixed 
point of the system (II, 2|) . The fixed point is not in H 5 ^ G and is given 
by 

(b) The system linearized about the fixed point has exponentially de- 
caying eigenvalues. 

(c) Every regular solution approaches the fixed point in I 2 norm. 

(d) Every solution blows up in finite time in H 5 / e norm. 

In a more technical companion paper [5] we study the solutions of (|1.2p 
after the time of blow-up in H 5 ^ 6 . Among the properties we prove that 
the H s norms for s < 5/6 are locally square integrable in time. Moreover, 
we study the long time behavior of (II. 2ft and prove that the fixed point is a 
global attractor. This is a consequence of anomalous or turbulent dissipation 
that was conjectured by Onsager [18] and related to Kolmogorov's prediction 
that in fully developed turbulent flow the energy spectrum in the inertial 
range is given by a power law 

(1.3) £(|£;|)= C oe 2 / 3 |£r 5/3 - 

In Section [6] of the present paper we observe that the energy spectrum 
of the fixed point of the forced model (|1.2j) reproduces Kolmogorov's law 
(II. 3p . Furthermore, regular solutions, i.e., solutions with bounded H 5 ^ 6 
norm, satisfy the energy equality, whereas the fixed point is not regular and 
does not satisfy energy equality. Moreover, anomalous dissipation cuts in 
for all the other solutions as well when the H 5 /® norm becomes unbounded. 

Organization of the paper. In section 2 we present a derivation of the 
dyadic model based on Littlewood-Paley operators. In section 3 we review 
existence and finite time blow-up results for the model obtained in the past 
few years. In sections 4 and 5 we prove the results stated above. In section 
6 we discuss how the dyadic model with the 3-dimensional scaling satisfies 
Onsager's conjecture and Kolmogorov's 5/3 law. 

Acknowledgements. The authors would like to thank Marie Farge, Jonathan 
Mattingly, Kai Schneider and Eric Vanden-Eijnden for very helpful discus- 
sions. S.F. was partially supported by NSF grant number DMS 0503768. 
N.P. was partially supported by NSF grant number DMS 0304594. 

2. Model based on Littlewood-Paley decomposition 

Let us start by considering 3D Euler equations with zero force (for the 
sake of simplicity): 

fin 

(2.1) — + (u-V)u + Vp = 0. 
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2.1. Euler equations on M 3 . Given a velocity u(x,t), let 



/ e~ ix< u(x,t) dx, £ G 
Jr 3 



be its Fourier transform. We shall use a standard Littlewood Paley decom- 
position. More precisely, we consider Fourier multipliers Pj (on _L 2 (R 3 )) 

so that their symbols Pj(£) are smooth and supported in |2 3 ' < |£| < 3(2 J ') 

and such that Pj(£) = Po(2~ J £) an d = Also we denote by the 

following multiplier Pj = Y%=-2 Pj+k- 

We apply the Leray projection T onto the divergence free vectors followed 
by an application of Pj on the equation (|2,ip to obtain: 

BP i] 

(2.2) Z£L +P .(T(u-Vu)) = 0. 

Our goal is to see what we can say about ||Pju||^2 because that is the relevant 
quantity in the definition of the Sobolev norm 

(2.3) \\f\\ H s = [ £ (1 + 2^)11^/111/ 

\j=-oo 

We use Bony's paraproduct formula for the nonlinear term. More precisely, 

Pj(T{u • Vu)) = N j>lh + N jM + N j>hh + N loc , 
where the low-high part is given by 

N i,lh= E Pj({Pj'u)-P 3 Vu), 
i'<i-4 

the high-low part is given by 

N JM= E P,((P,n)-P,vVn), 
i'<i-4 

the high-high part is given 

N jM = PjiiPjiu) ■ Pj'Vu) + Pj({Pj'u)-PjiVu), 

j'>j+4 j'>j+4 

and the local part is given via 

Ni oc = Yl Pj({Pfu) ■ P r Vu) + Pj({Pfu)-P r Vu). 

j-4<j'<3+4 i-4<i'<i+4 

Now let us start modeling. We recall that the Fourier multiplier Pjf(x) 
are given via 

pTf(o=pj(oho 

has the symbol Pj(£) which is, roughly speaking, supported for |£| ~ 2 J . 
This combined with the belief that only local frequency scales are relevant 



DYADIC MODEL OF TURBULENCE 



5 



in a turbulence cascade, motivates us to keep only local interactions in the 
nonlinear term, i.e. we keep only the equivalent of a modified version of N[ oc 
in the model that we propose below. More precisely, we introduce a model 
that describes the evolution of the coefficients 

(2.4) a J (t) = \\P j u(x,t)\\ L 2. 

Here Oj(t) represents the total energy in the shell. In order to introduce 
in the model an analogue of the Ni oc expression, we need to see what scal- 
ing comes out of \\PjU • V Pj/(u)\\ L 2 when j and j' are close. By Holder's 
inequality combined with the fact that Pj(£) are supported for |£| ~ 2 J we 
obtain: 



(2.5) \\PjU ■ \/Pjiu\\^2 < ||ij-tt||ioo HVij-'ttll^a 

(2.6) < 2^ j \\Pju\\ L 2 ■ 2 j '\\P j/U \\ L 2 

(2.7) = 'li'djdy. 

where to obtain (|2.6|) we use Bernstein's inequality which can be stated on 
M. n in terms of Littlewood-Paley operators as follows 

II-Fj^Hl 1 ? — ^*'~^ n: '\\Pju\\u>, for q > p. 

Now we are ready to propose the following model for the 3D Euler equa- 
tions: 

(2.8) ^ = 2*G-i) J_ 1 - 2§V;+i, 3 e Z. 
Note that due to fl2.4j) . the L 2 (R 3 ) norm of the fluid velocity u is 

1/2 



\\ u Wl 2 = 

The Sobolev norm (12.31) of u is 



E i 

i j=-oo 



(2.9) \\u\\ H > = £ (l + 2 2j ' 



2.2. Space periodic 3D Euler equations. Consider now Euler equations 
in a periodic box $7 = [0, L] 3 . By Galilean change of variables we can assume 
that the space average of u is zero. Given a velocity u(x,t), let 




= _ / e i ^r x - k u(x,t)dx, k G Z 3 , 



2tt , 

n 
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be its Fourier coefficients. Define a 2 (t) to be the total energy in the shell 
2> < \k\ < V +1 

(2.10) a j (t) 2 = J2 \Mt)\ 2 , J = 0,1,2,.... 

2J<|fc|<2J+ 1 

An analysis similar to that given in Subsection 12.11 motivates the following 
model of the 3D Euler equations: 

^ = 2*Ci-U J_ 1 - 2§ Vi+l, J = 1,2,... , 

(2.H) r 

— = -a a 1 . 

Due to (|2.10p . the L 2 (S7) 3 and H S (Q) 3 norms of the fluid velocity u can be 
defined as 

/ \ 1/2 / x 1/2 

/ oo \ / oo s 

(2.12) \\u\\ L 2 = \\a\\ P = [Y,\ a j\ 2 , \Hh° = ^2 2 ^|a/ 

\J'=0 / \j=0 

3. Global existence and blow-up 
We write the dyadic model (|2.1ip as 

^1 = A J ~ 1 a 2 _ 1 - X j aj a j+1 + fj, j > 0, 

(3.1) f f 

da 

— = -a ai + Jo, 
at 

with the forcing fj > for all j. 

For notational simplicity we adopt A- 7 as the scaling parameter in the 
computations in sections [3] and [H We do this to illustrate that the results 
are qualitively independent of the exact choice of lambda ( which depends on 
the spatial dimension and the construction of the model). As we discussed 
in sections 1 and 2 the relevant lambda for our 3 dimensional dyadic model 
is 2 5 / 2 . This exponent determines the values of the critical exponents in the 
Sobolev space results proved below. 

We say that a(t) = (ao(t), ai(t), . . . ) is a solution to (|3.ip if it is a classical 
solution in the usual sense. More precisely, we have the following. 

Definition 3.1. A solution on [0, T] (or [0,oo), if T = oo) of (J31J) is an 
Z 2 -valued function a(t) defined for t G [0, T], such that Oj £ C 1 ([0, T]) and 
a,j(t) satisfies (|3.ip on (0,T) for all j. 

Note that if a(t) is a solution on [0, T], then automatically a,j £ C°°([0, T]) 
for all j. The global existence of a classical solution for any initial data 
a(0) £ I 2 was proved in Cheskidov [3] for a closely related system. The idea 
of the proof is as follows, see Theorem 4.1 in [3] for more details. Given an 
arbitrary time interval [0, T], consider a sequence of Galerkin approxima- 
tions to (13. lh . It is easy to show that this sequence is weakly equicontinuos 
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on [0, T]. Therefore, thanks to Ascoli-Arzela theorem, we can pass to a 
convergent subsequence and obtain the existence of a weak solution to (|3.ip 
on [0,T]. Since the nonlinear term has a finite number of terms, this weak 
solution is also a classical solution. Hence we have the following result. 

Theorem 3.2. For every a £ I 2 and f £ I 2 , there exists a solution a(t) to 
(I3T1) on [0, oo) with a(0) = a . 



Define the H s norm of a as 

\a{t)\\ H s= [^2^af(t) X 

Due to (I2.12p . || • can be viewed as the model analogue of the Sobolev 
norm of a fluid velocity. For each scale j we define the energy as 

(3.2) Ej{t) = ^a){t), 

and the energy of the whole system 



oo 

(3.3) E (t) = -\\a{t)\\l=Y,E 3 {t). 

3=0 

We define the energy of a "box" Bj as 



(3.4) EB J (t) = ^E j (t). 

Formally multiplying (|3.ip by aj and summing over all j gives 



j=0 j=l j=0 j=0 

We recall that A = 2 5 / 2 . Hence, if ||a(i)llH 5 / 6 ^ s bounded on some time 
interval, then the summations on the right hand side of (|3.5p are uniformly 
convergent, and we obtain the following energy balance property: 

, CO 

(3.6) jMt)=Y.^r 

3=0 

This motivates the following. 

Definition 3.3. A solution a(t) of (|3.ip is called regular (or strong) on 
[Ti,T 2 ] if 1 1 a(t) 1 1 jys/e is bounded on [Ti,T 2 ]. 

Hence, the total energy of a regular solution to (|3.ip satisfies the energy 
equality (|3.6p . A regular solution also satisfies 

(3.7) j t E B{J) (t) = X J - 1 a 2 _ 1 aj, J > 0. 
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Treating the j-th equation in (j3. If) as an ODE for a,j(t) gives 



(3.8) 



where 



Oj(t) 

a (t) 



1 



1 

W) 



aj(0)+ / A^V^tMt) dr 



o 



J > 0, 



ao(0) + /o / Kr)dr 



o 



(3.9) fi(t) = exp ^ X j a j+ i(r) dr 

Hence if Oj(0) > for all j, we conclude that this property is preserved by 
evolution of (|3.ip and 

Qj(t) > for all j > 0. 

Thus, (|3.7|) implies that energy in each box .Bj of a regular solution increases 
monotonically in time, i.e. 

(3.10) ^E Bj (t)>0, J>0. 

The inequality (|3.1U|) implies that the system produces a successive cascade 
of energy into higher and higher scales J. 

The basic properties of system (|3.ip described above contribute to the 
proofs of certain results that have appeared in Katz-Pavlovic |12| . Friedlander- 
Pavlovic [10J, Kiselev-Zlatos [14] . We recall these theorems. 



Theorem 3.4. Let a G H s for some s > 5/2 with a° > for all j > 0. 

Then there exists a time T = T(||a°||#s) > such that a unique solution 
a(t) to HSU with a(0) = a exists and o(-) G C([0, T]; H s ). 



This theorem is proved in Friedlander-Pavlovic [lOj using Picard's fixed 
point argument in the case / = 0. In particular, this theorem implies local 
existence of regular solutions with initial data in H 5 / 2 . The local existence 
of regular solutions is not yet known for initial data not in H 5 ^ 2 . 

Theorem 3.5. Let a(t) be a solution of (I3.1I) with a(t) G H , aj(0) > for 
all j > ; and a(0) ^ 0. Then \\a(t)\\H s becomes infinite in finite time for 
all s > 5/6. 

The first proof of this type of blow-up result in H s for a system closely 
related to (13. ip with / = was given in Katz-Pavlovic [12]. Versions for H s 
blow up were given in Friedlander-Pavlovic |10| . Kiselev-Zlatoss |14] (the 
above theorem) and for a Navier-Stokes model in [12] and Cheskidov [3]. 

Crucial to these blow-up proofs is the cancellation of the infinite sums on 
the right hand side of (13. 5p leading to the energy properties (13. 6p and (|3.7p . 
In the context of the scaling of this present paper, the smallest value of s 
for which finite time blow-up is proved is s = 5/6. We note that the claim 
by Waleffe [20] to use such a proof to produce blow-up down to s > (for a 
variant of the model with 2 J scaling) is not justified because the telescoping 
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sums in (|3.5p do not converge for < s < sq, where so is a relevant critical 
Sobolev exponent (so = 1/3 for 2 J scaling). 

Consider now a forcing of the form /o > and = for all j > 1. 
Clearly, 

is a fixed point of (|3.ip . Since A = 2 5 / 2 , note that the fixed point is not 
regular, i.e., it is not in H 5 ^ 6 . Note also that it does not satisfy the energy 
balance equation (|3.6p . Indeed, for a fixed point we have 

d x 

Moreover, at the end of this paper we will prove that there is no global in 
time regular solution to (|3.1|) . i.e., every solution a{t) with a(0) G I 2 blows 
up in finite time in -£T 5 / 6 -norm. This result improves Theorem 13.51 

4. Uniqueness of the fixed point 
We continue studying the dyadic model 



da 4 

(4.1) 



—1 = X> -a] - \3 aj a j+1 + fj, j > 
da 

— = -a ai + /o- 

For the forcing of the form /q > and /j = for all j > 1 we will investigate 
the spectrum of the system linearized about the fixed point 

(4.2) M = {*-* + *VJo}- 
4.1. Uniqueness of the fixed point. 

Theorem 4.1. For every force f = (/o, 0, 0, . . . ), /o > ; there exists a 
unique fixed point of (|4.1j) 

(4.3) {a j } = {\-i+ 1 6^] r }. 
Proof. Consider 

(4.4) A j =\i-*f~*a j . 
Then the equations (|4.1|) for a fixed point become 

A)^ - AjA j+1 = 0, j> 0, 
A A, = 1. 

Clearly, 

(4.5) Aj = 1, j > 

is a solution to the above system, which corresponds to 

(4.6) a j = \-i + T> v r T . 
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Let us show that this is the only fixed point. Indeed, if {aj} is a fixed point, 
then 

(4.7) A A 1 = 1, ^ = £±i, j>l. 

From these equalities we obtain 

Ao = M _ A 2 A 3 A _ A 3 A 4 A _ A±A 5 

At V A 2 A Q A X ' A 3 AoAx' A a A Ai 

Thus, we have 

(4.8) ^. = Aj A j+1 , j>2. 
Therefore, 

(4.9) A j+1 = A Aj 2 , j > 2. 
Hence, 



(4.10) Aj 



4 +2j - 2)/3 A^- 2) , j>3odd, 
At 23 ~ 2)/3 Af~ 2 \ j> Seven. 



Finally, since A 2 = Aq, we obtain 

(4U) ^-{^-) + (l- 2 -)/3 ; .> 3even> 

Clearly, Aq has to be equal to 1 in order for {aj} to be bounded. Then 
Aj = 1 for all j. □ 

Note that the above argument can also be applied to a more general force 
of the form / = (/ , f u . . . , f k , 0, 0, . . . ) with f k > 0. 

Lemma 4.2. Let f = (/o, /i, . . . , 0, 0, . . . ) with /fc > 0. If there exists a 
fixed point {aj} of (14. ip . t/ien 

(4.12) aj = A5aj + i, yj >k. 

Proof. Let be a fixed point. Since > 0, we have that afca/u +1 > 0. 
Consider 

(4.13) ^ = J . 

Then = 1 and 

(4.14) A j A j+1 = A]_ 1 , j>k + l. 
Therefore, we have 

(4.15) ±t± = 4b±, j>k + l. 



DYADIC MODEL OF TURBULENCE 



11 




From these equalities we obtain 

A k __ A k+2 A k A k+2 A k+3 A k A k+3 A k+A A k A k+A A k+5 

A k+ i A k ' A k+2 A k A k+ i ' A k+3 A k A k+ i ' A k+A A k A k+1 

Hence, it follows that 

(4.16) ^ = AjA j+1 , j>k + 2. 



Aj+i = A k Aj 2 , j>k + 2. 

At^^ir', j>k + 3,j-k odd, 
A£-*- 2 WA{g>, j>k + 3,j-keven. 
A\, we obtain 

3-20- 2 )+(l+2^ 2 )/3 . . , . g . , 

A k , j > k + 3, j — k even, 

A k , j > k + 3, j — k odd. 

Clearly, A k has to be equal to 1 in order for {aj} to be bounded. Then 
Aj = 1 for all j > k, which implies that a* = \ l ^aj + \ for all j > k. □ 

Using this lemma one can prove the following theorem. 

Theorem 4.3. Let f = (/ , /i, ...,/&, 0, 0, ... ) with fj > for all j. Then 
there exists a unique fixed point {aj} of (|4. lfl . Moreover, 

(4.20) aj = \~iC, Vj > k, 
for some nonnegative constant C . 

4.2. Spectral stability. We consider the fixed point of (|3.6f) with / = 
(/o, 0, 0, . . . ), fo > 0. Rescaling the variables, we can assume that /o = A -1 / 3 
for simplicity. Then the fixed point is of the form {A _J//3 }. To investigate 
its spectral stability we write 

(4.21) a,(t) = A~3 +ebj. 
where the perturbation {bj(t)} satisfies 
(4.22) 

^ = \f(2\-hj-i - X~hj - b j+1 ) + e(A J '- 1 6|_ 1 - X j bjb j+1 ) j > 



(4.23) 



— = -\ sb Q -bi- eo Oi. 
at 
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The linearized equations for small e are 

(4.24) ^ = A^(2A-§6 i _i - X'hj - b j+1 ) j > 

(4.25) ^ = -X-ho-h. 

We seek a solution to (|4.24p - (|4.25p of the form 

(4.26) bj = Cj e^ 
with (jl real. 

Lemma 4.4. There are no positive eigenvalues to the system (|4.24p - (|4.25p . 
Moreover, there exist solutions in H s , s < 5/6 to the system (|4.24p - (|4.25p 
of the form (|4.26p with the eigenvalues [i negative. 

Proof Substitution of K26\\ into KM\ - K25\ \ gives 

(4.27) \~^c j+1 + ajcj - 2A~ic,_i = 0, j > 

(4.28) a c + A~^ci = 0, 
where 

(4.29) aj = A~^(A"4 + A - ^ /i), j > 0. 

We shall construct a sequence {cj}, Cj ^ 0, j > that solves the system 

(14371) - (pEgp . 

Define 

(4.30) = A - s2±l, j > 0. 

c j 

It follows from (|4T27]) and (j430|) that 
Hence 

(4.32) d n = [a n+ x,a n+ 2, ...], 

where [a n+ i,a n+ 2, ■■■] denotes the continued fraction 

K+l) «n+2, •••] = 



2 3/2 

23/2 a „ +1 + 



2 3/2 

23/2 a „ +2 + 



Observe that from (|4.29[) 



(4.33) lim a, = X'K 

jr'^OO 
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for any finite fi. Clearly for sufficiently large n the coefficients a n +m, m> 1 
in (|4.32j) are positive and the infinite continued fraction (14. 32ft is convergent 
to a positive number. Furthermore it follows that 



(4.34) doo = lim d n = A 2 < 1. 

n— >oo 

From (|4.27p - (|4.28|) we obtain the characteristic equation 

(4.35) - a = [ai,a 2 , ...], 

where ay's are given by (14.29p . We observe that (|4.35p has no roots fi > 0, 
since in this case all ay including ao are positive. 

Suppose there exists fi real and negative such that (|4.35j) is satisfied. 
Choose 

c = 1 
ci = X^d 

(4.36) 

j_ 

Cj = A« dj-idj-2---do. 

Then this sequence satisfies (I4.27p - (l4.28p by construction and in view of 
(Oil) 

(4.37) Cj for j » 1. 

Thus the I 2 norm of {cj} is finite. In fact the H s norm of {cj} is finite for 
s < 5/6. 

We are now going to study the characteristic equation (|4.35|) . Let X([i) 
denote the function 

(4.38) X(n) = a (fj,) + [«i(/x), a 2 {n), ■■■], 

with fi £ (—oo,0). Then to solve (|4.35p we look for zeros of the function 
X(ji). 

Since ctj(^) is given by (I4.29P we can write X(fi) in the form 
(4-39) X(,) = W(A4 + ,) + 23 ^ ) , 

From (|4.39p and (|4.29j) we observe that 

(4.40) X(0) > 0, 
and 

(4.41) lim XU) = -00. 

IX—*— 00 

It follows from (I4.40p and H4.41H that X((i) cannot be continuous and positive 
for all fi G (—00, 0). For fi G ( — A -1 / 3 , 0) each entry of the continuous fraction 
defining X(fi) in (I4.35P is positive. Hence X(n) is continuous in (— A -1 / 3 , 0). 
Let (/xojO) be the maximal interval of continuity of X{[i). Then X(p) is 
continuous at A -2 / 3 ^o- This implies that X(A~ 2 / 3 /xq) = 0. Figure 1 gives 
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Figure 1. Graph of X(n) 



a plot using Maple to graph the function X(n) and show that the zeros of 
X(n) are approximately \x = —1.4 and fi = —3.6. □ 



5. Behaviour of regular solutions 
In this section we show that every regular solution (i.e. with bounded 

#5/6 

norm) approaches the fixed point in the I norm linearly in time. 
Furthermore, we prove that every solution with initial datum in I 2 blows up 
in finite time in the i? 5//6 norm. 

Theorem 5.1. Let a(t) be a regular solution of (|4.ip on [Tx,!^] with f = 
(A -1 / 3 , 0,0,...). We write aj(t) = X~i +bj(t). Then 

||6(r 2 )|| P -||6(r 1 )||p<-i.A-*(r 2 -r 1 ). 

Proof. Note that b(t) satisfies the following system of equations: 



(5.1) 
dbj_ 

dt 
db 

dt 



AT(2A-%_! - X'Hj - b j+1 ) + (A^- 1 ^.! - X J bjb J+1 ), j > 1, 



-A 3 6 — b\ — 



DYADIC MODEL OF TURBULENCE 



15 



Multiplying it by bj and taking a sum from j = to j = k we obtain 



( 5 - 2 ) ^E 6 I = - A " l6 o-Mi-^i 



2dt 



j=0 



+ (A^(2A Vjj \bj - X'h) - bjbj+i) + (A' '/.j ./>,• - X j bp j+ i[ 



It now follows that 
i=o 



+ J2 X ^ ( 2A * b i-i&j " A ~^i - b j b j+i) - X k b\b k+1 . 



Also we can rewrite (|5.3[) as 
1 d 



3=0 



j=0 



j=0 



2X1 b k b k+1 - X k b 2 k b k+1 . 



However, 

A: 



J2x^{x~hj - xh j+1 ) 2 =J2x^(x~h 2 j + xh 2 j+1 - 2bjb j+1 ) 



3=0 



3=0 

-2 



k k 

1 ^ — "\ 2j_ q r — "\ 2j_ 

-X 3 2_^Xz bj + 2_^X3 bjbj + i 

3=0 3=0 

1 . r, , 2k | 1 



A i6g + A 3 + 3b\ +1 . 



Hence (|5.4p gives 



< 5 - 5 > !.tE'< 



2 

2 dt j 

j=0 



^A^"(A s&j - A^6 i+1 ) 2 + A~^6| 

j=0 



+ _AT+f - 2A^6 fc 6 fe+1 - X k b 2 k b k+1 . 
Now note that since a(t) is regular on [Ti,T2], we have that 

\ Si 



(5.6) 



3^coJ Tl 



lim / A^a-(t)dt = 0. 
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Hence, since bj = aj — X~^ 3 , we obtain 

(5.7) lim C 2 (h% + h 2 k+1 - 2\fb k b k+1 - X k b 2 k b k+1 ) dt 

= (T 2 - Ti) ( -A" - 2A"3 + X's 
= - l --\~ 1 HT 2 -T 1 ). 

Now ([531) and <[5T71) imply that 

(5.8) ||6(r 2 )|| P - ||6(ri)|| P < -\ ■ \~Ht2 - Tx). 



□ 



5.1. Blow-up in finite time. Using the results of the previous subsection, 
we prove that there is no global in time regular solution to (|4.1j) with / = 

;a 1 :i .o.o....). 

Theorem 5.2. All the solution to (|4.ip with f = (A _1 / 3 ,0,0, . . . ) blow up 
in finite time in H 5 / e -norm. More precisely, for every solution a(t) with 
a(0) £ I 2 , there exists a time 

0<t*< 2A3||6(0)|| Z2 , 

such that 

limsup ||a(i)||#5/6 = oo. 
t->t* 

Proof. First note that every ou(i) is continuous, i.e., a(t) is weakly in I 2 
continuous. Therefore, if ||a(t*)||fl-5/6 = oo, then 

(5.9) limsup || a(i) || #5/6 = limsup ||a(i)| | #5/6 =00. 

t-+t*- t^t* + 

Now assume that the statement of the theorem is not true. Then there ex- 
ists a solution a(t) to flU) which is regular on [0,T] with T = 2A 1 / 3 ||6(0)|| /2 . 
Then Theorem 15.11 implies that 

(510) ||6(T)|b 2 <||6(0)|b 2 -2- 1 A-lT 

= 0. 

Therefore, 

(5.11) a i (T) = A-i, 

i.e., ||a(T)||#5/6 = 00. Due to the remark in the beginning of the proof, 

(5.12) limsup ||a(t)||#5/6 = 00, 

a contradiction. □ 
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Note that there is no global in time regular solution to (|4.ip with a more 
general force / = (/o, 0, 0, . . . ), fo > 0. This follows by application of 
Theorem 15.21 to a rescaled system (14. If) . 

6. Onsager's conjecture and Kolmogorov's 5/3rd law. 

In the past 60 years since its original presentation much attention has been 
given to the statistical theories of turbulence developed by Kolmogorov and 
Onsager. For recent reviews see, for example Eyink and Sreenivasan [9], 
Robert [19]. Onsager [18] conjectured a clear distinction between 2 and 3 
dimensions: in 2 dimensions the energy of a turbulent flow is conserved, how- 
ever, in 3 dimensions dissipation of energy persists in the limit of vanishing 
viscosity. This phenomenon is now referred to as turbulent or anomalous 
dissipation. It is suggested that an appropriate mathematical description of 
3-dimensional turbulent flow is given by weak solutions of the Euler equa- 
tions which are not regular enough to conserve energy. Kolmogorov's theory 
[15j predicts that in a fully developed turbulent flow the energy spectrum 
in the inertial range is given by 

(6.1) E(\k\)=c e 2 / 3 \k\-^, 

where e is the average of the energy dissipation rate. This law stated in 
physical space means that the Holder exponent h of the velocity is 1/3 
in a statistically averaged sense. Onsager conjectured that for exponents 
h > 1/3 the energy is conserved and that this ceases to be true for h < 1/3. 
An elegant proof of the conservation of energy of weak solutions of the 3 
dimensional Euler equations in Besov spaces B^'°°, h > 1/3, was given by 
Constantin et al [6]. 

The model system that we study in this present paper reproduces the phe- 
nomenon described above. Regular solutions, i.e., solutions with bounded 
norni) satisfy the energy equality, whereas the fixed point 

(6.2) {a,} = {2-^+5/12^}. 

does not. Note that a necessary condition for a function in physical space 
to have Holder exponent h is that the spectral exponent is greater or equal 
to 1 + 2h (i.e., the spectral exponent 5/3 of the fixed point in our model is 
exactly the Onsager critical exponent corresponding to h = 1/3). 

The model (|1.2j) is derived under the assumption that a|(i) is the total 
energy in the frequency space shell 2i < \k\ < 2i +1 . Hence, by @2}, the 
energy spectrum for the fixed point is given by the expression 

(6.3) £(|A;|) = 2 5 / 6 /o|fcr 5/3 . 

Since the dissipation rate for the fixed point is equal to the energy input 
rate, we have 

e = a /o = 2 5 / 12 /o 3/2 - 
Thus the fixed point satisfies Kolmogorov 's law dEIJl with c = 2 5 / 9 . 
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In Section [5] we proved that every regular solution approaches the fixed 
point and, furthermore, blows up in finite time in H 5 / 6 norm. After this time 
(as we show in [5]) the solutions are not regular enough to conserve energy 
and anomalous dissipation for an inviscid system produces Kolmogorov's 
turbulent energy spectrum. 

Anomalous dissipation and "life after blow-up" has been examined re- 
cently in a linear 1-dimensional model where the system is "simple" enough 
to be fully solvable but still produces anomalous dissipation via energy cas- 
cade towards higher wave numbers [16]. In our more technical companion 
paper [5] we examine "life after blow-up" for the system (II. 2p . In particu- 
lar, we prove that the fixed point is a global attractor. Since the support 
of any time-average measure belongs to the global attractor, the average 
dissipation rate is equal to the dissipation rate of the fixed point. Thus 
Kolmogorov's law stated above for the fixed point is valid for the system. 

Remark 6.1. It is well known that the 1-dimensional inviscid Burgers equa- 
tion produces energy dissipation for weak solutions with shocks (as we 
noted, a rescaled version of (|1.2p can be motivated by Burgers equation 
in Fourier space). In contrast, the Cauchy problem for the 3-dimensional 
incompressible Euler equations remains open (and a major challenge). The 
3-dimensional Euler equations may (or may not) share with the model sys- 
tem (|1.2p the turbulent processes conjectured by Onsager for 3-dimensional 
fluids. The model per se could be written with n-dimensional scaling and 
the system itself has no behavior that is inherently 3-dimensional or in- 
compressible. We note, however, that the reduction of an incompressible 
3-dimensional vector model for the fluid equations to the scalar system (II. 2p 
given in [10] cannot be implemented in two dimensions. Furthermore, any 
such model that does not preserve enstrophy cannot be a suitable model for 
the 2-dimensional Euler equations. 
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